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QQ (57) Abstract: A method, computer program product and apparatus for estimating the pneumatic pressure of a tire on a vehicle, 
comprising the steps of: receiving as an input a vehicle status sensor signal; estimating first parameter values of an adaptive filter 
hased on a predetermined model on the vehicle status dependent on the vehicle status sensor signal; calculating, dependent on said 
first model parameter values, a first tire pressure indication value (6) being dependent on and indicative of the pneumatic pressure of 
the tire. Specifically, the vehicle status signal may be an input a wheel angular velocity signal (1) and a lateral movement indication 
signal; the adaptive filter is based on a predetermined model on lateral vehicle dynamics dependent on the angular velocity of said 
wheel and the lateral movement of said vehicle. 
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Tire Pressure Estimation 



Field of the invention 

The present invention relates generally to a system for tire pressure estimation in a 
5 vehicle carried on tire fitted wheels, and more specifically to the detection and estimation 
of physical parameters dependent on the tire pressure. 



Backgroimd 

In the design of cars and other vehicles a great deal of attention is today directed to 
10 safety aspects. One such aspect of great importance is related to the tire pressure in the 
wheels due to the fact that a sudden flat tire can easily cause an accident. In addition, in 
areas with a great risk of being robbed at the road an involuntary stop in cormection with a 
change of wheels entails an increased personal risk for the drivers and the passengers. 
Another aspect more related to driving economy is based on the fact that a tire havmg a 
15 low air pressure is worn out faster than a correctly filled tire. An indication of the tire air 
pressure is therefore a means for avoiding unnecessary costs. 

_/ 

The state of the art 

Direct measurement of tire pressure exists in prior art but requires a comparatively 
20 expensive equipment integrated with the tire. Traditionally, methods indirectly measuring 
the tire pressure are based either on a frequency analysis or on a tire radii analysis. The 
frequency ansdysis may be performed either by computation of the frequency spectrum 
using e.g. a Fast Fourier Transformation (FFT) or by estimation of the model parameters in 
a linear model of the tire. Usually, the tire is modelled as a damped spring and the 
25 estimated model parameters are related to the resonance frequency, which in its turn 
depends on the tire pressure. 

A tire pressure system based on wheel radii analysis is disclosed for example in the 
US patent no. 5,922,949 to Nalcajima. Since the angular velocity iacreases when the wheel 
radius decreases, i.e. when the tire pressure decreases, the relative tire pressure is estimated 
30 by comparing the angular velocity with the vehicle speed. However, the wheel radius is 
also affected by wear and tear and depends on cruising states, such as cornering and 
braking. Wheel radii analysis generally works well in good conditions but is less accurate 
in case of slippage of the wheels. An estimate value of each individual tire pressure itself is 
not possible to obtain through this prior art technology since only the difference between 
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the wheel radii of different wheels are investigated. Disadvantages with methods based on 

■ 

wheel radii analysis are for exana.ple a lack of robustness with regard to different driving 
styles and loss of functionaiity during low friction or in split |i situations, i.e. when there is 
different friction for different tires. A further disadvantage is that the whole algorithm fails 
5 in case of malfiinction in one of the wheel speed sensors. Moreover, simultaneous pressure 
falls in the tires may in certain cases remain undetected due to the comparative nature of 
the wheel radii analysis methods. 

Tire pressure systems based on frequency analysis are disclosed for example in the 
patent documents no. EP 0 700 798, EP 0 925 960, arid U.S. 5,753,809. | 

10 In EP 0 700 798, the frequency analysis is performed using a Fast Fourier 

Transformation (FFT) on the extracted high frequency components. The result is a 
resonance frequeiicy that is compared to a predefined reference frequency. 

In EP 0 925 960, the wheel speed signal is filtered using a band-pass filter to 
achieve a wheel speed signal corresponding to the first or the second resonance peak ^ 

15 dependent on the wheel speed. In this piece of prior art a second-order discrete time model 
for the frequency analysis is used instead of the Fast Fourier Transformation. The discrete 
time model includes two unknown parameters that are identified by the least squares 
method. The resonance frequency is fiirther computed from the above mentioned two 
parameters and may fiarther be related to tire pressure. 

20 The U.S. patent U.S. 5,753,809 to Ogusu et al discloses a method estimating the 

tire pressure using a second-order linear prediction model for the rotational speed of the 
vehicle wheel and vibration of the tire using parameters identified from values of the 
rotational speed of each individually treated wheel. A resonance frequency is foimd from 
the identified parameters, and the tire pressure is estimated from a linear relationship 

25 between the pressure and the resonance frequency. This method decreases the 
computational complexity in comparison with the method using the Fast Fourier 
Transform, since less arithmetical processing is required. Low-capacity memory is 
satisfactory and the method may be implemented in a low cost system. This piece of prior 
art is directed to a tire pneumatic pressure estimating system and discloses the use of a 

30 wheel speed sensor positioned in each wheel. The sensor signals are gathered and 
combined in a signal processor, which outputs a tire pressure indicator signal for each 
wheel. This system is confined to producing a tire pressure signal within the accuracy 
limits of the wheel speed sensors. 

Disadvantages foimd in prior art tire pressure estimation methods based on 
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frequency analysis is a relatively weak correlation to pressure and a sensitivity to various 
sources of vibration e.g. dependent on the vehicle velocity. 

The US Patent No. 5,629,478 to Nalcajima et al discloses a method and device for 
detecting a pneumatic pressure drop of vehicle tires. This prior art is specifically concerned 
5 with the problem of correcting for influence on the pressure drop detection of different 
interfering factors such as speed, turning radius, firont/rear acceleration and lateral 
acceleration. The correction is carried out by substituting, in the calculation of the pressure 
drop, a velocity value, a firont/rear acceleration value, a lateral acceleration value and a 
turning radius value of the vehicle with a polynomial expression to determine a correction 
10 value. The polynomial expression is specifically selected such that an influence of the 

turning radius is exerted on all terms. This pressure drop evaluation method thus makes use 
of a plurality of tire pressure parameters that are all dependent on wheel radii analysis, or 
expressed differently, a rolling radius analysis. However, the compensating properties 
seem to be limited. 

15 In short, well perfomiing tire pressure estimation methods should compensate for a 

number of different phenomena such as comering, acceleration, traction torque and 
fiiction, which may cause a virtual or experienced wheel radius difference. Prior art 
methods all seem to be insufficient in handling high velocities, i.e. situations when the 
traction torque is large. 

It should be noted that a tire generally refers to thick piece of rubber or similar 
material fitted onto the wheels of wheeled vehicles. The mbber is usually filled with a gas, 
most commonly compressed air. This text applies to all kinds of tires, although the 
examples mostly refer to £ur filled rubber tires. 

Object of the Invention 

The object of the invention is to solve the problem of providing an improved 
method and system for estimating the tire pressure in the wheels of a vehicle. 
Aspects of the problem are; 

- to achieve a tire pressmre indication with an improved degree of accuracy; 

- to provide a robustness with regard to for example various driving styles, cornering, 
low fiiction driving cases and split jx situations; 

- to provide a detectability of simultaneous drop of pressure in the tires. 
A further aspect of the problem is to provide an improved wheel speed signal from 

existing wheel speed sensors. 
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Summar y of the Inventioii 

The object of the mvention is achieved by a method and a system comprising the 
deteraiining or estimation of parameter values in an adaptive filter based on a 

5 predetermined or pre-selected tire pressure model, thus resulting in an adaptively 

determined equation system for calculating a tire pressure indication value. According to 
an aspect of the invention, the model is preferably based on latereil vehicle dynamics 
dependent on angvilar velocity of the wheels of a vehicle and parameters indicative of a 
lateral movement of the vehicle. The estimation of model parameter values is calculated 

10 dependent on the input of parameter signals indicative of lateral vehicle dynamics 

preferably in the shape of an angular velocity signal fi:om one or a pluraUty of wheels and 
preferably a yaw rate signal, possibly in combination with or as an alternative a lateral 
acceleration signal., A tire pressure indication value, for example in the shape of a tire 
pressure offset value, which is dependent on and indicative of the pneumatic pressure of a 

15 tire is then calculated dependent on the model parameter values. In contrast with prior art, 
the invention estimates parameter values of a predetermined or pre-selected tire pressure 
calculation model such that the model is adapted to the current specific situation. The 
inventor heis realised that by estimating tire pressure model parameter values in a recursive 
filtering process, a large number of error sources are taken into consideration and their 

20 influence on the parameter estimate value, and in the next stage the error infiucence on the 
tire pressure indication value, is suppressed. 

In a further aspect of the invention, a first tire pressure indication value which is 
based on a first dynamical model is used in conjxmction with a second tire pressm^e 
indication value based on a second, different dynamical model for tire pressure estimation. 

25 This estimation technique taken in conjunction with a specifically designed evaluation 
process employed in an advantageous embodiment of the invention entails a further 
improved application of the the pressure indication. 

The parameters of the adaptive model is preferably calculated by means of a 
recursive filter, preferably in the shape of an recursive least, square filter or often even 

30 more advantageously a Kalman filter operating on a state space model. 

An important aspect of the invention is embodied in complex models and filters 
integrating and comprising parameters representing dififerent basis in terms of physical 
principles or dynamics, such as combinations of lateral and longitudinal dynamics or 
addition of other kinds of parameters such as motor parameters and road friction. Different 
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sensor signals are integrated in these complex models thus embodying a sensor fusion. 

A specific embodiment of the invention applies a combination of first estimation 
model in the shape of an improved frequency analysis and a second estimation model in 
the shape of an improved wheel radii analysis based on measiued wheel speed signals. 
5 According to an advantageous aspect of the invention, the frequency analysis comprises a 
second-order discrete time model provided with a Kalman filter for spring parameter 
estimation, a feature that efficiently attenuates noise and increases the theoretical as well as 
the practical accuracy of the analysis. More specifically, this embodiment comprises a 
calculation of the friction parameter slip offset constituting a tire pressure indication value. 
10 The slip offset parameter is independent of friction and traction force, and therefore the 
invention provides a tire pressure estimation that is insensitive to different fiiction levels 
and is robust to dijfferent driving styles. An embodiment of the wheel radii analysis 
comprises a computation of: 

- Firstly, the relative difference in wheel radii of the front and rear wheels for the left and 

15 the right wheel pair, respectively, which gives an indicator of whether either one of the left 
or the right wheels has an abnormal pressure. Preferably, the relative differences in wheel 
radii are computed from a fiiction model designed to compensate for different phenomena 
such as cornering, acceleration, traction torque and friction, which may cause a virtual or 
experienced wheel radius difference. 

20 - Secondly, the relative difference in wheel radii of the left and right wheels for the front 
and rear wheel pair, respectively, which similarly gives an indicator of whether either one 
of the front or the rear wheels has an abnormal pressxire. This second wheel radii analysis 
is based on a sensor fusion model that in addition to tlie signals from the wheel speed 
sensors takes a signal from a yaw rate gyro. 

25 Embodiments of the invention may also comprise a wheel radius parity check in 

order to provide an indicator of whether one of the tires has an abnormal radius. 

According to further aspects, embodiments of the invention comprises a pre- 
processing provided with a calibration algorithm based on an auto-regression model for the 
measured wheel pulse signal in order to generate an improved, computed wheel speed 

30 signal for the tire pressure estimation. In particular, this is an advantageous featui*e in 
embodiments using a wheel speed sensor that comprises a toothed wheel and an edge 
detector coupled to each respective wheel. Sensors of this kind are common in the standard 
anti-lock bralcing systems (ABS) of modem vehicles. A typical problem with such a 
toothed wheel is that neither the distance between the teeth nor the size of the teeth ai e 
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imiform and consequently an inaccurate wheel speed signal is generated. This problem 
may get even worse during operation due to the fact that the toothed wheel generally is 

■ 

badly protected from mechanical damage, for instance from stones or gravel. Embodiments 
of the invention therefore comprise a correction method for correction of the inaccuracy 
5 emerging from the toothed wheel. 

An advantage with the present invention is that the functionality works well in 
vehicles having any number of wheels, even in one- wheeled vehicles and perhaps more 
practically in two-wheeled vehicles such as motorcycles, in standard four-wheeled cars or 
lorries with extra wheel pairs. 
10 Embodiments of the present invention are advantageously applied for wheel 

imbalance detection, surface texture estimation and surface unevemiess detection. These 
applications are frirther e^lained below. 

The invention is applicable for and adaptable to all kinds of vehicles carried on 
pneumatic tires, such as cars, motorcycles, lorries or trucks having any number of wheels 
15 in accordance with the different embodiments and further aspects described in the detailed 
description. 

Brief Description of the Drawings 

The present invention will be further explained by means of exemplifying 
20 embodiments in conjunction with the accompanying drawings, in which: 

Fig 1 shows a block diagram illustrating the general method steps as well as the 
functional means of an embodiment of the invention; 

Fig 2 shows a block diagram illustrating an embodiment of the invention applied in 
a four-wheeled vehicle; 
25 Fig 3 shows a block diagram of an embodiment of the sensor signal pre-processing 

shown in Fig 1 ; 

Fig 4 shows a toothed wheel of a wheel speed sensors used in connection with tlie 
invention; 

Fig 5 A shows a diagram of the angular wheel velocity signal output from an ideal 
30 toothed wheel; 

Fig 5B shows an exemplifying diagram of the angular wheel velocity signal output 
from an unideal toothed wheel; 

Fig 6 shows a diagram of an energy spectrum of an imcorrected angular wheel 
speed signal; 
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Fig 7 shows a diagram of offsets in radians plotted for the tooth of a toothed wheel; 
Fig 8 shows a diagram of the resulting energy spectrum of an angular wheel speed 
signal after correction in accordance with the invention; 

Fig 9 shows a general block diagram of the toothed wheel correction in accordance 
5 with the invention; 

Fig 1 0 shows a diagram of a exemplifying energy spectrum in the frequency 
domain for an angular velocity signal of a car; 

Fig 1 1 shows a flow chart of the separation of frequency bands in an angular 
velocity signal; 

10 Fig 12 shows a flow chart of the inventive reduction of the sample frequency in the 

lower frequency band; 

Fig 13 shows diagrams of frequency spectra of signals from different steps of the 
reduction process in the lower frequency band; 

Fig 14 shows a flow chart of the inventive reduction of the sample frequency in the 
15 higher frequency band; 

Fig 15 shows diagrams of frequency spectra of signals from different steps of tlie 
reduction process in the higher frequency band; 

Fig 1 6 shows a general flow chart of a data decimation functionality using aliasing; 
Fig 1 7A shows the input and output signals of the data decimation using aliasing; 
20 Fig 1 7 B illustrates interpolation; 

17 C shows a flow chart of an interpolation process; 

Fig 1 SCshows a diagram of resonance frequency estimation using a Kalman filter 
for three different tire air pressures; 

Fig 1 8 A and 1 8B show flow charts of embodiments of the adaptive model 
25 estimation; 

Fig 1 9 illxistrates an embodiment with an adaptive filter estimation based on a 
combination of yaw rate analysis and wheel radii analysis; 

Fig 20 shows the geometrical relations in four wheeled vehicle; 
Fig 2 1 and 22 show diagrams of resulting tire pressvire estimations with and without 
30 a non-linear filter for speeding up adaptation; 

Fig 23 shows the relations between offset values; 

Fig 24 and 25 show vector diagrams for detecting a low pressure tire; 

Fig 26 and 27 show embodiments for high precision yaw rate and road friction 

based tire pressure estimation; 
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Fig 28 — 32 show diagrams relating to a motorcycle application of the invention; 
Fig 33 and 34 show etnhodiments for integrating high precision yaw rate and 
vibration analysis in tire pressure estimation. 

5 Detailed Description of Embodiments 
General Setting * 

The invention compriseis a method and a reaUsation of the method that may be 
realised in hardware, software or a combination thereof. The most feasible realisation of 

> 

the invention is likely to be in the shape of a computer program product preferably 
10 comprising a data.carrier provided with program code or other means devised to control or 
direct a data processing apparatus to perform the method steps and functions in accordance 
with the description. A data processing apparatus running the inventive method typically 
includes a central processing imit, data storage means and an I/O-interface for signals or 
parameter values. The invention may 2ilso be realised as specifically designed hardware 
15 and software in an apparatus or a system comprising mechanisms and functional stages or 
other means carrying out the method steps and functions in accordance with the 
description. In practice, the invention will typically be implemented in a digital or 
analogue sensor control system of a vehicle. Such a system includes wired, optical or 
wireless data communication links for communicating the output of a sensor to a control 
20 unit. The control unit itself is provided with a data processor of the above mentioned kind. 

» 

The input and output of different functional blocks is in this text described in terms of 
signals, by which should be understood analogue signals, sampled signals or data 
representing the described entities or parameters. 

25 Tire Pressure Estimation 

The tire, pressure estimation is calculated based at least partly on the angvdar 
velocity of each wheel subject to estimation. The angular velocity signal is preferably 
taken from angular velocity sensors of an anti-lock bralce system (ABS) installed in the 
vehicle or from sensors specifically installed for the purpose of the invention- Such 

30 sensors, however, are generally not very accurate for example due to imperfections in the 
manufacturing and to mechanical wear. An embodiment of the invention is therefore 
provided with a pre-processing stage comprising a calibration of the measured sensor 
signal, a feature which greatly enhances the accuracy of the computed wheel speed signal. 
Embodiments of the inventive pre-processing and calibration algorithms, respectively, will 
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also be further explained below. In cases where a sufficiently accurate angular velocity 
signal is provided, the invention is embodied without the pre-processing stage. 

In one embodiment a first tire pressure estimate is calculated based on a first 
model of vibration phenomena, and a vibration analysis is carried out using an adaptively 
5 estimated model of the firequency spectiaun of the vibrations. A first tire pressm-e indication 
value is then calculated dependent on parameter values resulting firom the vibration 
analysis. In order to further improve the accviracy of the tire pressure indication value, the 
firequency based tire pressure estimation is combined with a second tire pressure estimate 
based on a second model, preferably involving an analysis of the wheel radii also based on 

10 the angular velocity signal. 

Fig 1 shows a flowchart illustrating the functional components as well as method 
steps of an embodiment of the invention. An angular velocity signal 1 or cd is received 
from each wheel in a pre-processing means 2 carrying out for example a calibration 
algorithm. Embodiments of the pre-processing also comprises a data decimation algorithm, 

15 which decreases the computational burden with a factor of up to 5 compared to a non-pre- 
processed signal or data sequence. A pre-processed and thereby improved angular velocity 
signal 3 is then transmitted to a model estimation means 4, which is devised to estimate 
parameters 4 upon which the tire pressure depends based on said first physical principle of 
vibrations. Preferably, a model estimation is carried out based on an individual firequency 

20 analysis of each wheel. A first tire pressxire estimate 7 is then calculated firom the 
parameter values 4 obtained in step 4. 

In parallel, the original 1 or the pre-processed 3, as shown with the intermittent line, 
angular velocity signal co is also transmitted to a second tire pressure estimation means 8, 9 
for tire pressure estimation based on said second physical principle. Preferably, the second 

25 tire estimation comprises a wheel radii analysis means 8 devised to calculate the wheel 
radii and a second tire pressure calculation means 1 5 producing a second tire pressure 
estimate 16. In addition, different relations or estimates 9 of the respective wheel radii are 
possibly calculated to be used as a complement to the first tire pressure estimate 7 of step 
6. Preferably, the relative difference in wheel radii of the firont and rear wheels for the left 

30 and the right wheel pair, respectively, is computed firom a friction model. The relative 
difference is used as an indication of abnormal pressure in any of the right or left wheels. 
In a further developed embodiment of the invention another input signal to the wheel radii 
analysis is taken firom a yaw rate gyro and combined with a similar computation of the 
relative difference in wheel radii of the left and right wheels for the front and rear wheel 
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pair respectively. The combination with the yaw rate gyio signal is used to compensate for 
cornering, acceleration and Mction, and gives as a result an indicator of whether any of the 
front or the rear wheel has an abnormal pressure. Possibly, a wheel radii parity check may 
also be included in order to generate an indicator of whether one of the four tires has an 
5 abnormal radius. 

The wheel radii analysis and embodiments using Kahnan filtering in the 
model estimation (explained below) is preferably based on a firiction model, wherein the 
slip s of a wheel is defined as the relative difference of a driven wheel's circumferential 
velocity, cowrw, and its absolute velocity Vwthus s = (cDwr^- Vw) / Vw In these equations a>w 
10 is the wheel angular velocity and Tw the wheel radius. The absolute velocity of a driven 
wheel is computed firom the velocity of the two non-driven wheels using geometrical 
relations. The firiction coefficient p., is defined as the ratio between the traction force, Fi, 
and the normal force, N, on the driven wheel. Said firiction coefficient is also referred to as 

■r . 

the normalised traction force. The slip slope, k, is defined as d\x/ds for |li=0, and is different 

15 for different surface textures. Said slip slope is commonly referred to as the longitudinal 
stiffiiess, but in this text slip slope is preferred. However the slip s is not zero when tHe 
traction force, Fi, is zero, and due to a small difference in effective wheel radii a slip offset, 
5, appears. This slip offset, 5, has the advantageous featm e of being independent of firiction 
and traction force, and as is only dependent of the tire radii. The wheel angular velocity co 

20 is computed based on the signal firom the wheel speed sensor and the slip s is in its turn 
calculated from the angular velocity. Furthermore, the engine torque M required in some 
embodiments of the invention for the computation of a normalised traction force [x, is 
preferably calculated based upon signals firom engine sensors. For the model estimation , 
explained below, the parameters co and p., are then input into an adaptive filter, e.g. a 

25 Kalman filter, adaptively estimating the slip slope, k, and tlie slip offset, 5. 

Returning to Fig 1, the first tire pressure estimate 7 as well as the second tire 
pressure estimates 1 6 and possibly also the wheel radii computations 9 are then received in 
an evaluation means 1 0 devised to evaluate the input signal and to calculate a resulting tire 
pressure estimate 1 1 . Thereby, first tire pressure estimates 7 from driving situations which 

30 are handled with uncertainty in the firequency based model estimation are corrected or 
compensated for by means of the wheel radii dependencies on the actual tire pressure. 
Examples of such driving situations compensated for are comering, acceleration and 
firiction. Described in a different wording, the evaluation means 1 0 is in fact a means for 
fusion of different sensor signals or different signals derived from the same sensor in order 
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to malce an optimised use of the available infoimation. Experimental tests have shown that 
a tire pressure estimation deviating as little as about 15% from the real tire pressure is 
achievable with the invention. 

As finally shown in Fig 1, the resulting tire pressure estimate 1 1 is according 
5 to one exemplifying embodiment transmitted to a change detection means 12 devised to 
generate an alarm signal in response to detected changes according to predetermined rules. 
For example, a first alarm signal may be generated if the tire pressure of one tire decreases 
more than 25% in a short time interval indicating a puncture, and a second alaim may be 
generated if the tire pressure increases above a predetermined level, 

10 

Four- Wheeled Vehicle Application 

In Fig 2, an embodiment of the invention is shown by means of a block diagram 
representing method steps as well as functional units of the invention. The shown 
embodiment is designed for tire pressure estimation in a four wheel vehicle, and it should 

15 be appreciated that embodiments adapted for other numbers of wheels would have the 
same basic structure. This embodiment comprises or utilises a wheel speed sensor 
1 0 1 , 1 02, 1 03 , 1 04 detecting the angular velocity for each wheel. The wheel speed sensors 
are preferably existing equipment in the vehicle, for example provided in connection with 
an anti-lock breaking system (ABS). 

20 The signal from each wheel speed sensor is individually pre-processed by sensor 

signal pre-processing means 111,112,113,114 and is then filtered by an adaptive or 
recursive filtering means 121,122,123,124 adapted for a frequency model estimation and 
calculation of parameter values upon which the tire pressure depends. A firequency analysi; 
is then carried out in frequency analysis means 13 1,132,133,134 on the obtained parametei 

25 values and generates for each wheel a calculated tire pressure signal, which is transmitted 
to a first logical unit 140. In the first logical unit 140, the calculated tire pressure signals 
are processed e.g. with regard to predetermined threshold values and tire pressure signals 
for the firont left wheel FL, the right left wheel RL, the front right wheel FR and the rear 
right wheel RR, respectively, are generated and output from the first logical unit 140. 

30 Possibly in addition or altematively, one or several waining signals Wl is generated and 
output from said fust logical unit 140. 

A fixrther development of the invention may also comprise a second logical unit 
142, into which the tire pressure signals FL,RL,FR,RR are input. The second logical unit 
142 is optionally also coupled to a means for slip offset estimation 144 devised to perform 
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a slip offset estimation based on the original angular wheel velocity signal from wheel 
speed sensors 101,102,103,104 or the pre-processed angular velocity from the pre- 
processing means 111,112,113,114. The input of angular wheel velocity is in Fig 2 for the 
sake of simplicity shown as co. The estimated slip offset values for each wheel are 

5 transmitted to the second logical unit 142 as an input to a further improved evaluation and 
estimation of the tire pressures. Similarly, there may also be provided wheel radii analysis 
means 146 or wheel radius offset estimation means 148 adding information for improved 
estimation of the tire pressures. The second logical iinit 142 is adapted to generate a second 
set of improved tire pressvire signals 150 for the wheels or a warning signal W2 dependent 

10 on predetermined rules. 

Pre-processing 

The purpose of the pre-processing stage is to enhance the quality of the signal upon 
which the model estimation is based. An embodiment of the pre-processing means adapted 

15 for this specific type of measurement device is described as an example. Other types of 
sensors appUed in the invention would need a similar general pre-processing scheme, but 
the details would be adapted to the specific correction problems in each case. Fig 3 shows 
an embodiment of a pre-processing method as well as the functional means for carrying out 
tiie method. An angular velocity signal 202 from a wheel speed sensor is thus input ta a 

20 pre-processing means 201 which comprises a sequence of a detector signal correction- 
means 204, a re-sampling means 206, a band pass filtering means 208 and a down- 
sampling/data decimation means 210. 

Detector Signal Correction 

25 The most common angular wheel velocity measvirement device in modern standard 

cars and other vehicles is based on a toothed wheel arranged on and rotating wifli the 
vehicle wheel and a sensor, e.g. an optical or an inductive sensor, detecting each tootli 
edge. The time elapsed between the detection of the two last tooth edges is measured and 
an angular velocity is calculated with the equation: 

30 (/.) fi) = Q.n IN)-tst [rad/s] 

where N is the number of teeth on the toothed wheel and At is the elapsed time between 
two detected teeth. A problem with this type of devise is that the toothed wheel is usually 
not ideal for example due to mechanical tolerances in the production, i.e. the distance 
between the teeth and the sizes of the teeth are not completely uniform and thereby the 
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angle increment between all teeth is typically not unifonn. Furthermore, the toothed wheel 
is often insufficiently protected against mechanical damage and deformation for example 
from stones or gravel. These and other phenomena make the angular velocity signal 
incorrect. 

5 Fig 4 illustrates a toothed wheel 40 1 wherein the outline of an ideal wheel is drawn 

with an intermittent line and a more realistic, not ideal wheel is drawn with a fiill line. The 
a,- is the real angle between teeth i and i-1 and di is the offset error. Fig 5 A shows the 
angular velocity, as a function of the time t and for a car travelling with a constant velocity, 
the theoretical output signal © from a wheel speed sensor haviag an ideal toothed wheel 

10 and thus resulting in a constant co. Fig 5B shows similarly an example of an output signal 
CO from a real, imperfect toothed wheel sensor, resulting in a periodically constant © witti. 
the periodicity N, i.e. the number of teeth. In the frequency domain, the periodic 
disturbances in the angular velocity signal cause an impulse with large number of 
harmonic components. Fig 6 shows an energy spectrum of an uncorrected angular wheel 

15 speed signal, in which spectrum the components originating fi'om the imperfect toothed 
wheel dominate and makes other frequency components unobservable. The interesting 
parts of the signal often have the same frequencies as the disturbances. 

The period of the periodically constant velocity signal is the same as for the vehicle 
wheel itself, i.e. with a normal wheel having a circumference of approximately 1,8 m and a 

20 vehicle velocity of 80 km/h the wheel rotates on lap in T/ap = 1.8/(80/3 .6) = 0.08 Is. This is 
equivalent with a frequency of flap = 1/0.081 =12.3 Hz. This should be compared with the 
theoretical fact that the imperfections from the toothed wheel have a fundamental 
fi-equency at about v/(3.6*L8) Hz, thus depending on the actual velocity v. As plotted in 
Fig 6 the fundamental frequency of the impulses is at about 1 1 .6 Hz, which corresponds 

25 well to the tiieoretical frequency of the toothed wheel. Fig 6 thus shows an the frequency 
spectrum, i.e. the Foimer transform of an uncorrected wheel speed signal and illustrates the 
distortion caused by mechanical errors. 

According to the invention, the imperfections of the toothed wheel are identified and 
thereafter the disturbances originating from the toothed wheel are eliminated fi'om the 

30 angular velocity signal. Again referring to Fig 4, the offset errors 5, are identified in order 
to find the correction factors. The real angle a/ is unknown but can be approximately 
calculated with: 

(2) Of; = 07 • At 
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where co is the current angular velocity of the wheel. The current wheel speed ay can be 
calculated in different ways, and in order to keep the algorithm causal only old 
measurements are used in accordance with the invention. In one embodiment © is 
calculated using tiie average velocity during the last lap, according to 

5 (3a) fl> = 2;r/(/;-^/_;y). 

Preferably, an adaptive algorithm, e.g. a general Recursive Least Squares (RLS) or a 
Kalman filter algorithm is used to identify the offset error. The advantage with using an 
adaptive algorithm like RLS is that if new errors occur it adapts to the new toothed wheel. 
With this approach the angular velocity is expressed as: 
10 (3b) o>=(27t/N+©i)/At. 

The linear regression model for identifying the errors of the toothed wheel pertaining to 
mechanical causes is: 

0 = iS, 5^ - s^Y 

<p(fi ) = (0 1 • • • 0) position i mod N is\ 
(4) (o{t,) = 27tl(t^-t,_^) 

N 

z 

In equation 4 co is the current angular velocity, N is the number of teeth and t, is the time 
15 when the last tooth passed the sensor. The current angular velocity is not available and is 
therefore calculated from the past data. The past data can be taken firom an arbitrary 
number of revolutions (k), but there is a trade-off betweep sensitivity and delay for the 
signal. As an additioned constraint the sum of the whole offset must be 0.- 
The solution of a linear regression model with RLS is: 

y(t) = <p^ (J)9 + e(^t) regression model 



20 (5) 



Kit) = FitMO 



->] 



A+<p' (t)p{t-l)<p(t) 

The offset errors (5i 62 • • • 6n) are then estimated by combining equations (4) and (5) and 
applied them to the signal from the angular velocity measurement unit. 

TjTjical offset errors estimated by means of the described algorithm are shown in 
Fig 7, the offset error in radians is plotted for each tooth. As seen in Fig 7, some teeth have 
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relatively large offsets which distort the wheel speed signals quite prominently. When the 
offset errors are identified the wheel speed signal is corrected by for example: 

N 



(6a) 



^corr 



wherein the angle is corrected with the identified mechanical errors before the correct 
wheel speed signal is calculated. 

In a second alternative the wheel speed signal is corrected by: 
(6b) 27r 



At- 



wherein the time elapsed between the last two teeth is corrected with the knowledge of the 
mechanical errors before the correct wheel speed signal is calculated. 

10 In the equations 6a and 6b here At is the original signal, 5 is the offset for the actual 

teeth and cd current wheel speed using calculated e.g. with equation 3. 

In Fig 8 the resulting energy spectrum, i.e. the Fourier transform of the wheel speed 
signal after correction is presented, and it is clear that almost all harmonics due to 
distortion from mechanical errors are totally eliminated and the real signal components can 

15 be extracted. 

Fig 9A shows a flow diagram summing up the described embodiment used in a 
variety of the detector signal correction 204 of Fig 3, thus taking as an input the elapsed 
time between two teeth At upon which a linear regression model is created 902. Based on 
the linear regression model the toothed wheel is identified with an RLS and the offset 

20 errors are estimated 904. Thereafter, a correction of the wheel speed signal is carried out 
906 and a corrected elapsed time Atcorr is output. Similarly, Fig 9B shows the alternative 
correction process used in a alternative variety of the detector signal correction 204 of Fig 
3, inputting At into a step of creating a linear regression model 908, identifying and 
estimating the offset due to mechanical error 910, correcting the wheel speed signal and 

25 outputting a corrected cocorr. 



Re-sampling and Reduction of Sampling Frequency 

The re-sampling process 206, band pass filtering 208 and the down-sampling 210 in 
the pre-processing as described in Fig 3 is described in conjianction with Fig 10 to 17. 
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The most feasible application of the invention is in a digital control system where 
the angular velocity signal (o(t) from a wheel speed sensor is sampled and the sampled 

signal or data ©(nTi), with the sampling frequency of — , is available to the system 

components. As explained above, this is sufficient information for a tire pressure 
5 estimation to be independently estimated for each wheel, in particular in a low pressure 
interval. When carrying out the frequency analysis there are basically two pressure 
dependent frequency intervals from which information about the tire air pressure can be 
extracted. These two intervals are here referred to as the low frequency interval and the 
high frequency interval, respectively. Fig 10 shows an exemplifying energy spectrima in 
10 the frequency domain for an angular velocity signal of a car, wherein the low frequency 
interval 1001 and the high frequency interval 1002 are indicated. 

The low frequency component is substantially due to the fact that the side of the 
wheel is elastic, and consequently a vibration may occur if the corresponding frequencies 
are excited. This frequency is highly dependent on the tire tj^e, since different tire types 
15 have different elasticity. The high frequency substantially originates in the tire tread, and 
since the size of the tire tread is much smaller than the wheel, the fi^equency from the tire 
tread is high. 

In one embodiment of the invention, information from both intervals are used in 
order to reduce the sample frequency and thereby also reduce the data rate. Therefore, the 

20 frequency information from the two frequency intervals is separated by means of two band 
pass filters. Fig 1 1 shows a flow chart of such a separation wherein a sampled angular 
velocity signed co(nTi) is input into a first band pass filter 1101 having a low pass band in 
the interval [fi f2] Hz and into a second band pass filter 1 102 having a high pass band in 
the interval [fs f4] Hz. The output from the lower band pass filter 1101 is a first new signal 

25 coi(nTi) and from the higher band pass filter is a second new signal 0i)2(nTi), Thereafter, th< 
lower frequency interval spectrum is moved to [0 fa-fi] Hz and the higher frequency 
interval spectrum is moved to [0 fi-fa] Hz in, and in connection with the frequency move a 
reduction of the sample frequency is carried out. 

Fig 12 shows a flow chart of the steps for moving the frequency spectrum and for 

30 reducing the sample frequency in the lower frequency band. The sampled signal G)i(nTi), 
shown in Fig 13 a, is input into a data decimation stage 1201 which decimates the data with 
a factor ki by saving every ki :th sample and discai'ding the rest of the samples. The 
decimation results in a signal yi(nTiki), which is plotted in Fig 13b. The signal yi(nTiki) i: 
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thereafter modulated in stage 1202 with cos(27iflTl ki) thereby generating the modulated 
signal y2(nTiki), which is plotted in Fig 13c. In order to enable further decimation the 
signal YzCnTiki) is filtered through a low pass (LP) filter 1204, thereby avoiding aUasing. 
The result of the low pass filtering is a signal in y3(nT|ki) plotted in Fig 13d. Finally, the 
5 signal ysCnTiki) is in stage 1206 decimated by a factor ki, by saving every k2:th sample anc 
discarding the rest of the samples. The thus output signal firom this data decimation process 

has a sample firequency that is reduced from — in the input signal to — ^ — of the resulting 

output signal y4(nTikik2) plotted in Fig 13e is thereafter available for processing. 

Fig 14 shows a similar flow chart of the steps for reduction of the sample frequency 

10 in the higher frequency band. The sampled signal (02(nTi), shown in Fig 15a, is modulated 
in stage 1401 with cos(27cf3Tl) thus using the frequency f3 and thereby generating the 
modulated signal y5(nTi), which is shown in Fig 15b- The signal y5(nTi) is then filtered foj 
example through a low pass (LP) filter 1402, removing the component which will cause an 
alias for the decimation step. The result of the low pass filtering is a signal in y6(nTi), 

15 plotted in Fig I4c, which is input into a data decimation stage 1403 which decimates the 
data with a factor ka by saving every ky,th sample and discarding the rest of the samples. 
The decimation results in a signal y7(nTik3), which is plotted in Fig 15d. The thus output 

signal from this data decimation process has a sample frequency that is reduced from— in 

the input signal to — ^ of the resulting output signal y7(nTik3) and is available for 
20 processing. 

The modulation teclmique described above in connection with the data reduction 
schemes is computationally rather expensive, in particular for the high frequency interval 
where a great deal of signal processing has to be carried out at the original sampling 
frequency. An alternative embodiment of the invention is therefore based on a utilisation o 

25 the aliasing phenomenon. This embodiment is shown by way of an example applied on the 
high frequency interval. Fig 1 6 shows a flow chart wherein the high frequency angular 
velocity signal CD2(nTi) is roput in a data decimation stage 1601, which decimates the data 
with a factor k4 and outputs the resulting signal y8(nTik4) by saving every k4:th sample an( 
discarding the rest of the samples.. The input signal o>2(nTi) as well as the output signal 

30 y8(nTilc4) are shown in Fig 17a and 17b, respectively. 
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Interpolation of Non-equidistant Data 

When signals from for example the ABS in a pre-processing stage have been corrected 
regarding imperfectiohs in the toohted wheels the data is interpolated in order to enable use 
of standard theory for equidistantly sampled systems. The idea behind interpolation is to 
5 convert data iD^ (ti ), 02 if 2 \ ca^ (^). . . collected at time , ^2/3 • • • to a new data sequence 

Q>[ Qq \ a>2 (fo + t\ 0)2 {tQ t 2T), . . . This process must be performed online and several methods 

such as linear interpolation or spline interpolation are possible. Fig 17B illustrates the idea 
behind interpolation, i.e. comprising the steps to: 
1 . Find the last sample taken at time < nT 
10 2. Find the first sample taken at time > nT 

3. Calculate the gradient between said last and said first values 

4. Calculate the value at time nT using the interpolation along the straight line 

The implementation is straightforward and one altemative for linear interpolation is 
illustrated as a flow chart in Fig 17C, comprising the steps of 1702-1719. 

15 

BP-filtering to prepare signal for model identification 

In embodiments identifying a low order parameter model, such as AR(2), to estimate 
the tire pressure in a vibration based model, the resonance frequency is preferably 
separated from other disturbances. This is achieved using a digital band pass (BP-) filter 
20 expressed in the form 

(7) y{n) = ai y (« - 1) + . . . + a y{n - A^)+ b^x^n) + 61 jc(n - 1) + . . . + jc(rt - Af ) 

where ^(«) is the original signal and y(n) is the filtered signal. The characteristics of the 

filter are decided by the filter coefficients 

(8) ajfe ^ = l...^' 

25 and 

(9) bt / = 1...M 

There are many different choices of the filter coefficients and they are calculated 
using per se known standard methods. Examples are discrete time Cauer, Butterworth, 
Chebyshev I or Chebyshev II filters among others. In order to reduce the filter order a 
30 Cauer filter may however be prefered compared to Butterworth, an example of 
considerations to take account of. 

Model Estimation 
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The vibration analysis of the invention is based on a spring model of the tire 
and on the well known fact that road roughness as well as engine torque excite resonance 
frequencies in the tire. Using wheel velocity signals from wheel velocity sensors as a 
starting point, a distinct peak in the energy spectrum dependent on the tire pressure is 
5 foimd for example at the frequency interval 40-50 Hz. An estimated tire pressure is 

* 

therefore calculated by means of a frequency analysis on wheel velocitj' sif^^l*^ ^ 
known manner. 

The tire is modelled as a damped spring, which expressed in a discrete time 

form is: 

10 (10) ^r=-^l^^-l--^2^/-2+^/- 

With the model applied in a tire air pressure system yt is the actual angular velocity of the 
wheel and et is external noise from the road and the present state of the vehicle. 
This is a second order auto regression model with a transfer function 

1 



(11) G(q)^ 



1,, — 



15 where q is an operator such that q" yt = yt-i. The purpose of the model estimation, i.e, the 
system identification, is to estimate the model parameters ai and ai from which the tire 
indication can be calculated. There are alternative methods for carrying out the estimation 
of the parameters, however, it is important that the algorithms are executed on-line, 
meaning that every sample should be treated immediately rather than as a batch of data 
20 collected and processed off-line in order to avoid time delay in the sensor system. 
According to the invention, the estimation is preferably carried out by means of an 
adaptive estimation algorithm such as a Recursive Least Squares method (RLS) or a 
Kalman filtering. 



25 Recursive Least Squares Model Estimation 

In an embodiment of the inventive model estimation, the model parameters are 
estimated by means of a recursive least squares method (RLS). This embodiments 
comprises the introduction of the parameter vector and regression vector as: 

<9 = (a, a^f 



(12) 



30 to simplify the notation. Then the criterion 

(13) 
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is miniinised with respect to dy by means of the RLS. The solution to (9) is obtained by 
means of the recursive algorithm: 



(14) K(t)^P{t)q>{t) 

^ = - 1) + K(tj^{t) - <p ^ (/) h - 1) 

The design variable in RLS is the forgetting factor A A large X (close to 1) makes the 
5 algorithm to remember a large number of old measurements and therefore has a low 
adaptation gain. A small X ,on the other hand, makes the algorithm faster, but it is instead 
more sensitive to disturbances. Both parameters ai and a2 have the same adaptation gain, 
and only one of them is useful for the tire pressure estimation. Experimental tests have 
shown that the parameter ai is strongly correlated with the pressure, while a2 has a poor 
10 correlation. 

Fig 1 8 A shows a general flow chart of the inventive model estimation using a 
Recursive Least Square, comprising the following steps: 
1 802 establishing or creating a linear regression model; 

1 804 estimating parameter values by means of a Recursive Least Square process; 

15 1812 transforming parameters to a resonance frequency. 



State Space Model Estimation 

Another embodiment of the model estimation and parameter estimation is 
based on Kalman filtering based on a state space model. In the inventive tire air pressure 
20 calculation system, the state space model is for example expressed as: 

(15) x,^{a^ a:,Y 



A, =1 



10 <ii) 

The design variable in the Kalman filter is the covariance matrix O/, which is a 
description of the relation between Qt and J?/. The Kalman filter allows for different 
adaption gains {qi and qi^ for the parameters al and a2, respectively. , which is an 
25 important advantage in the model estimation since one of the parameters takes into account 
the different tire air pressures and the other parameter handles some of the external 
disturbances. 
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After the estimation of the parameters values, the parameters are transformed or 
calculated into a resonance frequency. With the transfer function expressed in equation 
(11) above, the tire model of equation (10) can be written as: 

5 The poles in G(q) is a complex conjugate and the angle a between the poles and the real 
numbers axis correspond to the resonance frequency. The calculation of the angle a is in 
an exemplifying embodiment in the tire pressure calculation system implemented with 
computational advantages as: 



(17) Qj = arctan(- 

10 Since a is in the range [0 ic] it has to be transformed with respect to the sampling period Ts 
in order to obtain the resonance frequency oOres, where 

(18) fl,^=-^[Hz]. 

When the tire air pressure decreases, the dependent parameter ai decreases and as a 
consequence also the resonance frequency ©res decreases. 
15 Fig 18B shows a general flow chart of the inventive model estimation using a 

Kalman filter, comprising the following steps: 
1808 establishing or creating a state space model; 

1810 estimating parameter values by means of a Kalman filtering process; 

1812 transforming parameters to a resonance frequency. 

20 Fig 1 8 shows a diagram of the course of the resonance frequency estimation with the 

Kalman filter for three different tire air pressures, viz curve 1=2,3 bar, curve 2=1,95 bar 
and curve 3=1,6 bar. 



Tire Pressure Indicator Using High Precision Yaw Rate 

25 In one embodiment of tiie invention an adaptive Tire Pressure Indicator (TPI) is 

achieved by combining a first tire pressure estimation based on a first physical principle, 
preferably being a wheel radii analysis on rolling radius phenomena, with an estimation 
based on a second physical principle, preferably being an analysis on yaw rate phenomena. 
The wheel radii analysis is preferably computed dependent on signals from wheel angular 

30 velocity sensors of an ABS sensor, and the yaw rate analysis is preferably computed 
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dependent on signals from a gyro. In an advantageous variety of this embodiment, signals 
from a lateral acceierometer is additionally used to enhance the performance of the system. 

Fig 1 9 shov\rs a flov^chart schematically illustrating the functional 
components as well as the method steps of an embodiment of a tire pressure indicator 
5 comprising an adaptive model estimation 1902 receiving as input signals a yaw rate signal 

\j/ from a gyro and wheel angular speed signals ©i from an ABS. The model estimation 
1902 is preferably realised as an adaptive filter, for example a Kalman filter based on a 
space state model £is described above. The input signals are possibly pre-processed in a 
pre-pi*ocessing stage 1901 in order to improve signal quality and usability. The input 

10 signals are also received in an observability evaluation stage 1906 devised to determine the 
observability of the parameters and to transmit an enable signal to the adaptive model 
estimation stage in case the observability is positively determined. The adaptive model 
estimation stage 1902 produces as an output offset signals 5front and Jrear, respectively, and 
Sfioni hat and Srear hah respectively. The offset signals Sfroi^t and Srear are input into a non- 

15 linear filter 1904 which in a feedback connection transmits its output g/ront and qrear back to 
the adaptive model estimation stage 1902 in order to speed up the adaptation £ifter a 
pressure drop. The q —parameters describe the speed of the adaption for the two offset 
values and is normally small (slow). However, if a pressure drop would occur the 
adaptation speed is increased in order to enable tracking of the new situation. The values oi 

20 signals representing S/ronthat and Srear hat are transmitted to an evaluation stage 1908 devised 
to evaluate tire pressxire changes and to generate change indicator signals, for example in 
the shape of an alarm in accordance with predetermined rules. 

Fig 20 shows a schematic drawing of a four-wheeled vehicle during comering, in 
which drawing the geometric relations are defined for the wheel velocities used to compute 

25 a yaw rate estimation from wheel speed signals. The wheels are in this application denoted 
rl for rear left, rr for rear right, ^ for front left and fr for front right. A coordinate system 
indicating the x, y, and z-directions is also drawn in Fig 20. 

In the adaptive model estimation 1902 a high precision yaw rate value is calculated 
by means of a first yaw rate signal from the gyro combined with a second yaw rate value 

30 calculated dependent on wheel angular speed signals. The following equations are used in 
this process. 
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(19) ^ 
wherein is the yaw rate firom a gyro; 

Vx is the velocity of the vehicle in the x-direction taken from an ABS system; and 
Oy is the acceleration in the y-direction optionally taken from a lateral accelerometer, and 
5 input into the adaptive model estimation. 

The curve radius is defined as the distance between the rotational centre of the car 
and the middle of the rear axis, and is computed according to the following relation, where 
R is defined as the distance to the center of the rear wheel axle, 

V R ^ + ^ 

(20) = = / ^ 

10 The angular wheel velocities o) for each of the respective wheels are received from 

an ABS and the inverse R'^ of R is solved in order to avoid numerical problems in certain 
driving cases, e.g. driving straight ahead. This results in 

(21) R"^ — ^ ^fr _^ ^rr^rr 

where the wheel radius is denoted ri. 
15 The wheel radii ratio is subject to an offset: 

(22) ^-1 + ^^ 

r 
rr 

The influence of the offset on the denominator is negUgible, so in preferred 
embodiments the following expression is used for inverse curve radius: 



(23) R-' = 



1 2 



^(l + ^^o.)-!^ 



~ "I ^n,^ 



rr ^rr 



20 wherein the computable quantity 



(24) RZ' = 



1 2 f^_^^ 



is used for the inverse curve radius. 
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The velocity at the centre of the rear wheel axle is 
(25) v,^= '* ^ " r 

where r denotes the nominal wheel radius. 
Finally the yaw rate is computed with 

5 (26) 

In a practical implementation of this embodiment, the sensor measurements 
are typically as follows: 

- j^iCt) is from a yaw rate sensor, i.e. gyro signal; 

- Mt) = v,^„.^i?:;„,, , from rear ABS sensors; 

10 - ^wlnrar IS computed as K'^ above; 
V as V : 

- 3^3(t) = V:,^nufront > ^^^^ ^BS sensors; and 

- Rl]/ront and V, „ are computed in the same way as R^^ and v, ,„ , where is 
replaced with o)^ and co^ is replaced with o . 

15 AH these sensor measurements are subject to an offset and measurement noise given by 
the relations: 
(27) 

1 2 ^aiO 

1- 1 ^ 

where Srear is an offset that depends on the relative tire radius between rear left and right 
20 wheels and Sfront on the relative radius between the left and right front wheels. S/ront is 

defined in the same way as Srear- 

In accordance with the invention, the offsets are estimated adoptively by recursive 
least squares (RLS) algorithm, least mean square (LMS) or a Kalman filter. 
The Kalman filter is preferred and has the advantages over RLS in the following way: 
25 - Firstly, a model for the variation of the true yaw rate can be incorporated in the Kalman 
filter. 
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Secondly, different time-variations of the sensor offsets can be used. For instance, 

temperature can influence the gyro offset variation, a cold start can make the filter 

forget more of the gyro offset than the ABS offset 

The Kalman filter is completely specified by a state space equation of the form 

x{t'hl)=Ax{t)-hBv{t) 
j;(0=Cx(/)+6(0 

where the covariance matrices of v(t) and e(t) are denoted Q and R, respectively. 

Hie problem when designing an implementation is to set up the state space model 
An exemplifying embodiment uses the state vector: 



5 (28) 



(29) 



^(0 



gyro 



rear 
\^ front J 



and a continuous time state space model is: 
(30) 
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0 



1 0 0 
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rt 



x^t^rear 



L CO., , 0) 

" -I- 1 rr 



CO 



rr 



0 



0 
0 



1 



€0 



x,m, front 



fr 



J 



It is here assumed that there is an unknown input v(t) that affects the yaw acceleration, 
which is a common model for motion models, basically motivated by Newton's law F = 
15 ma. 
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A discrete time state space model is described ia the following equation (31) and is 
used by the Kalman Filter. 



(31) 
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K0= 
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0 
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1 2 



CO, 



x^m^rear 



0 
0 



CO 



rr 



1 0 0 



0 



1 



6) 



x^m^front 



if)^e{t) 



J 



The relative wheel radii offsets do not vary during normal driving with correct 
pressure and therefore the adaptation gains for the wheel radii offsets are low. To speed up 
the adaptation after a pressure drop in the rear wheels a non-linear filter is applied, in 
accordance with: 



(32) 



g(t)rear =^ ^ tf g(0«a, < 0 

if giO > h, alarm and increase one sampel 



SiOrear = 0 



Where^(0^,, =;^2(0- 



1 



CO 



rl 



CO, 



^4(0 



0) 



rr 



J 



In the same way a non-linear filter can be applied for a pressure drop in the front 
wheels. 

The system monitors 5^^^^ and 5j^^^^ if the any of these values exceeds a threshold h an 

15 alarm is provided. The following predetermined rules are used to decide which tire 
pressures that are incorrect. 
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Event 


Alann 




Pressure drop in rear nght wheel 




Pressure drop in rear left wheel 




Pressxire drop in front right wheel 




Pressure drop in front left wheel 



Fig 21 shows a plotted curve of the results of the tire pressure estimation 
without a non-linear filter, and Fig 22 shows the results with a non-linear filter for 
5 speeding up the adaptation gain. 

In a stepwise description of this embodiment for tire pressure indicator, the 
inventive method comprises the following steps (referring again to Fig 19)- 

1) Collecting measurements: 

(a) yaw rate signals from a yaw rate gyro; 
10 (b) angular wheel speed signals from an ABS Sensors; 

(c) possibly lateral acceleration signal ay from a lateral accelerometer; 

(d) possibly wheel steering angle from a wheel steering angle sensor; 
input to a pre-processing stage 1901 or into an adaptive model estimation step 1902. 

2) Pre-processing and filtering of raw sensor data, comprising: 
15 (a) Scaling of sensor inputs to physical constants, 

(b) Low-pass filtering of yaw rate to reduce quantization and noise error effects for 
instance averaging over a few samples. 

(c) Rotational synchronisation of cog stamps to avoid cog deformity error effects 
i.e. calculate wheel angular velocity by using a fixll revolution of the wheel. 

20 (d) Similar treatment of sensor signals Ic) and Id) to reduce known error effects for 

instance low-pass filtering. 

(e) Performing simple diagnosis on sensors to take care of sensor built-in diagnosis 
functions for instance sensors usiag a zero level voltage to indicate internal failure. 

3) Calculating filter inputs and parameters in error model, comprising: 

25 (a) Calculating the inverse curve radii estimate from front axle wheelage (equation 

24, vsdth wheel index changed to coji and co/r,) 

(b) Calculating the inverse curve radii estimate from rear axle wheelage (equation 
24). 
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(c) Calculating the vehicle velocity estimate from ABS sensors (equation 25). 

(d) Calculating the yaw rate estimate from front and rear wheelage (equation 

26). 

(e) Calculating the wheel radii offset error propagation function for rear axle. 
5 (f) Calculating the wheel radii offset error propagation ftinction for front axle. 

(4) Examining data quality by given norms to produce statistical matrices for Kalman 
Filtering, for instance: 

(a) Low velocity increases noise and other errors in yaw rate from ABS sensors 
estimated => increasing the covariance matrix R. 
10 (b) Standstill car assures yaw rate to be exactly zero => decreasing the 

covariance matrix R. 

(5) Applying the Kalman Filter equations. 

(6) Detecting changes using a non-linear filter. 

(7) Output adaptive filter 

15 (a) Providing relative wheel radii S^^ between left and right wheels on rear and 

front wheelage to be used for a tire pressure indicator systems. 
If S^Qs > ^ exceeds the threshold h a pressvire drop alarm is provided. 



Tire Pressure Indicator Using High Precision Yaw Rate and Road Friction Indicator 
20 In one embodiment the tire pressure indication is achieved by combining 

parameters based on lateral vehicle dynamics with parameters beised on parameters 
dependent on road friction in an integrated tire pressure model. Preferably, the lateral 
dynamics parameters are derived from a high precision yaw rate computation as described 
above and the road friction parameters is derived from a torque signal indicating the torque 
25 of the vehicle motor together with a gearing ratio indication signal. The integrated tire 

pressure model is in accordance with the invention realised in an adaptive filter or adaptive 
estimation algorithm such as a Recursive Least Squares Method or a Kaknan filter, as 
described above. 

Using the geometric vehicle relations as described in connection with Fig 20, a 
30 method for calculating the road friction dependent on angular wheel speed is described in 
the following equations. 



Symbol 


Description 




Angular velocity, driven wheel [rad/s] 
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Angular velocity, non-driven wheel [rad/sj 




Wneel radius, driven wheel [mj 


r„ 


Wheel radius, non-dnven wheel [mJ 




Engine torque [Nm] 


I 


Gearing ratio from engine to driven wheels [-] 


N 


Normal force at driven wheel [N] 


s 


Wheel slip [-] 


M 


Normalized traction force [-] 



10 



The longitudinal stiffness k , which during normal driving conditions can he modeled as 



(33) 



where 



TJ 



2r,N 



is the normalized traction force and 



(34) 



is the wheel slip. Here N is the tire normal force that depends on the mass, the vehicle 
geometry and the vehicle state (such as the current velocity and acceleration/retardation). 



In practice one must take into account that the wheel radii and r„ are unknown. By 
introducing 6 as the relative difference in wheel radii, i.e., 



(35) 



r r 



15 (36) 



and approximating the wheel slip as 



CO 



-1, 



n 



one obtains the model 



(37) 



1 



s » 



s^S ^—ju-^S 
k 
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The parameters in this model, Ijk and S , can be estimated from measured slip s^^ and 
traction force fi .using an adaptive filter (such as a recursive least sq\iares algorithm or a 
Kalman filter) and the state space model 

5 where = (l/A: SY ^ ^(ju, l/ , and and are process and measurement noise, 
respectively. 

This method is in one embodiment used separately for the left and the right side of 
the vehicle. This yields two relative radii offsets 



(38) 



10 (39) 



rl 



^ right 



fr 



rr 



assuming a front wheel driven car. In case of a rear driven car the fractions on the right 
hand side need two be inverted- The relative wheel radii offsets do not vary during normal 
driving with correct pressure and therefore the adaptation gains ( q^^j^ and q^^g^^ ) for the 

wheel radii offsets are low. To speed up the adaptation after a pressure drop in the rear 
1 5 wheels a non-linear filter is applied. 

giOi^t =Oif gif)uji < 0 

if g(t)i^ji > hy alarm and increase q^^j^ one sampel 

giOieft = 0 

where s(t)f^jt = y^^j^ - {/^tjeji^iOi + ^(Oi)- the same way the non-linear filter can be 

applied for the right side. The system monitors S^^j^ and S^^,^ if the any of these values 

exceeds a threshold h an alarm is provided. The following logistics are used to decide 
20 which tire pressures that are incorrect. 



(40) 



Event 


Alarm 




Pressure drop in front left wheel 


^ieji <-h 


Pressure drop in rear left wheel 


bright > h 


Pressure drop in front right wheel 
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Fig 23 shows how the relative wheel radius offsets from high precision yaw 
rate computation and road friction computation are related. Fig 24 and 25 show how these 
offsets can be used to find out if a wheel has a low tire pressure, in relation to the other 
5 wheels, and also to identify which wheel that deviates in tire pressure from the other ones. 
More specifically Fig 24 shows an illustration of how offset vectors in a normal tire 
pressure situation sxmi up to zero, but for a small estimation error due to the fact that the 
tires are not absolutely similar. Fig 25 shows a situation where the rear right tire has a low 

pressure and {J^^,^ and S^^^ are significantly larger than the estimation error. 

10 Fig 26 shows schematically the offset computation from CAN signals 2603 

inputting yaw rate and angular wheel velocity to a high precision yaw rate computation 
stage 2601, and torque and gear indication to a road friction indication computation stage 

2602.The output offset values 2604 from the two computation stages are dj^j^ and S^^i^^^ . 

Fig 27 shows schematically the offset values 2604, 2701 being input into an 
15 integrated tire pressure model as mentioned above comprising parameters representing the 
respective offiset values relating to lateral vehicle dynamics and road friction. In an 
integrated tire pressure model for example reaUsed by means of a Kahnan filter, the output 
2703 may be adjusted to indicate a pressure warning 2704, a the indication or 
identification 2706 and a confidence indication 2708. 

20 

Tire Pressure Estimation for Motorcvcles 

One application of the invention is directed to tire pressure esthnation in 
motorcycles. In different embodiments of the invention different sets of parameters and 
input sensor signals are combined or integrated in the tire pressure estimation model. Fig 
25 28 shows the geometric relations for normal and lateral forces on a motorcycle expressed 
in well established notational terms, and Fig 29 shows in the same manner the notation for 
angular velocities. 

In the motorcycle embodiment of the invention the angular velocity of the 
front and the rear wheel, for example detected by sensors in an ABS equipment, are input 
30 to a predetermined tire pressure estimation model, possibly in conjunction with additional 
sensor signals for improving the estimate. An additional parameter used in advantageous 
embodiments is the roll angle as detected by a roll angle sensor. The roll angle is in a 
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simpler embodiment used as a threshold to turn off indications j&om the tire pressure 
estimation system when the roll angle is outside a permitted interval. In another 
embodiment the roll angle is comprised and integrated in tlie model in order to increase the 
roll angle interval in which the tire pressure estimation is xxseful. 

5 Fig 30 shows schematically an embodiment of the invention taking as an 

input angular velocity signals for the front wheel 3001 and for rear wheel 3002 
respectively, and the roll angle 3003 to a tire pressure estimation stage 3008 comprising an 
adaptive model as described above in this description text. The angular velocity signals are 
before their input into the estimation stage pre-processed in a pre-processing stage 3006 

10 preferably comprising a combination of a first and a second low pass filter, a resampling 
and a data decimation stage as described in the pre-processing section of this text. 

In a more advanced embodiment, schematically shown in Fig 3 1 , further 
single or combined signals from sensors detecting parameters such as torque 3012, brake 
3014, air temperature 3016 and tire temperature 3018 information are comprised in an 

15 integrated tire estimation model 3020 to improve the system further. Similar to the 
previous embodiment, the additional sensor signals may be processed in a specifically 
adapted pre-processing stage3013, 3017, 3019 in order to improve signal quality. 
In addition to the above equations, the model adapted for motorcycles takes into account 
the following basic equations for detecting changes in effective wheel radius between the 

20 front and the rear wheel. The wheel slip for the rear wheel is: 



s = 



II, 



(OjVj I r 



CO. 



CD 



Sj. -1 



If no torque is applied to any of the wheels, the slip are equal to zero which gives 



CO. 



CD 



= 0=>^J 



CO 



CD, 



25 The performance is improved by removing outliers, which typically are connected to 
different types of driving such as acceleration, turning and braking. The parameter S,. is 
preferably estimated using a recursive filter such as recursive least square (RLS) or 
Kalman filter. 

An improved embodiment comprises an algorithm integrating torque 
30 information in the following manner. A measurement of the applied torque is used to 

calculate the applied force between the wheel and the ground. When the nomial forces on 
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the wheels are known, the normaUsed traction force // = — is estimated. The relation 

N 

between slip and normalised traction force is illustrated in the diagram of Fig 32. 
Preferably the system is designed for small slip values and thence the curve can be and is 
in one embodiment made linear according to 



5 (43) k = ^ 



US 



Because the true velocity of the rear wheel is approximated by the angular velocity 
multiplied with the wheel radius of the front wheel this curve has a significant offset 

(44) ^^ = ^1^.0 

The linear model for the slip slope is then modelled as 
10 ju = k{s-s),la order to estimate the parameters k and S , a good regression model required, 

preferably 
(45) 

A state space model for an adaptive filter such as a Kalman filter is expressed as: 

^ ^ yif)=Hm)^e(f) 
where 

Q{f)= Ev{ty if) 
R{})= Ee{f)e^ if) 

(47) KO'^^meamnK/ 

H{t)=m 1) 

.(o=(iA(/) 5{f)r 

In order to calculate the normalised traction force, the normal forces acting on the 
motorcycle are estimated. Derivation of these forces are straightforward, especially if the 

20 system takes account of the roll angle by being tumed off when leaning in tums, as in the 
first altemative mention above. The alternative method is to extend and integrate the roll 
angle in the model to make it able to handle leaning of the motorbike. An embodiment 
realising the latter altemative will have to handle the lateral forces acting on the wheels, 
which in their turn influences the longitudinal stiffiiess. In accordance with embodiments 

25 of the invention employing a Kalman filter, this is efficiently dealt vsdth in a properly timed 
Kalman filter since the longitudinal stiffiiess can change very fast e.g. due to different 
types of roads snow, ice etc, but the wheel radius offset is changing slowly. 
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In order to compensate for normal and lateral forces in turns models for the forces 
acting on Hxe wheels are derived in accordance with the following equations, for the 
notation referring to Fig 28 and 29. 

Force acting on front wheel normal to tihie road is expressed as: 

5 (48) Njr ^j^mg-——(uza •^l.uO) 

Force acting on rear wheel normal to the road is expressed as: 

(49) N^=-£-rr,^^—{^^-l^ue) 

Force acting on the rear wheel in the road perpendicular to the direction of travel is 
expressed as: 

10 (50) = .^^^^ +i^ue) 

Force acting on the front wheel in the road perpendicular to the direction of travel is 

expressed as: 

rc^\ „ m sin O / . , a\ 

(51) Fj^ ^^^^[^uzo ^Ifue) 

The inventive tire pressure estimation applied for motorcycles is thus adapted 
15 with the aid of the algorithms of this section. 

Fusion of Wheel Radii and Vibration Based Tire Pressure Estimation 

The embodiment of the invention is directed to tire pressure estimation by 
combining wheel radii based and vibration based methods in an integrated tire pressure 

20 model. More specifically, this embodiment employs the results of the different 
embodiments described above. 

Fig 33 illustrates schematically thcfimctional stages of this embodiment, taking 
CAN signals 3301, such as angular wheel velocity 3302, torque 3304, gear ratio 3306 and 
yaw rate 3308 as an input into a model integrating a high precision yaw rate estimation 

25 3310 with a road friction based estimation 3312. The output results in terms of a first 
pressure warning signal 3316, a first tire identification 3318 and a confidence indication 
3320 are taken as an input to an integration model in a fusion algorithm, such as a Kalman 
filter as described above. Furthermore, input signals in the shape of angular velocity 
signals 3326 are taken from e.g. ABS sensors 3324 are input inot a vibration based 

30 estimation model 3328 as described above. The vibration analysis output in the shape of a 
second pressure w^aming 3330, a second tire identification 3332 and a second confidence 
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indication 3334 are input into the integration model 3322. The integration model comprises 
an adaptive filter for estimating integrated model parameters and thereafter calculating a 
resulting pressure warning signal 3336 and a resulting tire identification 3338. 

Fig 34 shows a flow chart for the method of generating a warning in accordance 
5 with the steps of 3402 to 3416, using the signals for pressure warning Wi , tire ti and 

confidence c,- where i = 7,2. First there is a check if any of the algorithms signals an alarm. 
If yes, there is a check if it is the same tire. If yes, the sum of the confidence is compared to 
a predetermined threshold. If yes, an alarm is sent that tire tj has low pressure. If the 

methods have warnings for different tires ^ algorithm compares the confidence 
1 0 and the method with largest confidence is then checked to a predetermined threshold. If the 

confidence is larger than the threshold there is an alarm for that tire. 

The invention has been described by means of exemplifying embodiments and it 

should be appreciated that several designs and combinations are possible within the 

inventive concept as defined in the claims. Furthermore, all different aspects and functional 
15 features of the invention may appear independently applied. 
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Claims 

1 . A method for estimating the pneumatic pressxjre of a tire on a vehicle, comprising the 
steps of: 

5 receiving as an input a vehicle status sensor signal; 

estimating jBrst parameter values of an adaptive filter based on a predetermined 
model on the vehicle status dependent on the vehicle status sensor signal; 

calculating, dependent on said first model parameter values, a first tire pressme 
indication value (8) being dependent on and indicative of the pneumatic pressure of the 
10 tire. 

r 

2. The method as recited in claim 1, fixrther comprising the step of combining said 
adaptive filter model with a rolling radius based calculation of a second tire pressure 
indication value dependent on a vehicle status sensor signal representing wheel angiilar 

15 velocity, 

* 

3 . A method for estimating the pneumatic pressure of a tire on a vehicle as recited in 
claim 1 , comprising the steps of: 

receiving as an input a wheel angular velocity signal (1) and a lateral movement 
20 indication signal; 

estimating first parameter values of an adaptive filter based on a predetermined 
model on lateral vehicle dynamics dependent on the angular velocity of said wheel and 
the lateral movement of said vehicle; 

calculating, dependent on said first model parameter values, a first tire pressure 
25 indication value (5) being dependent on and indicative of the pneumatic pressure of the 
tire. 

4. The method as recited in claim 3, wherein said lateral movement indication signal is 
generated by a sensor indicatuag the yaw rate of said vehicle, said sensor preferably 

30 being a gyro. 

5 . The method as recited in any of the preceding claims, further comprising a lateral 
movement indication signal generated by: 

integrating, in an adaptive filter, parameters relating to lateral dynamics and 
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parameters relating to rolling radius phenomena; 

estimating the values of said integrated parameters dependent on a lateral 
movement indication signal generated by a sensor indicating the yaw rate, preferably a 
gyro, and on a wheel angular velocity signal; 
5 calculating an improved lateral movement indication signal dependent on said 

integrated parameter values. 

6. The method as recited in any of the preceding claim, further comprising the step of 

integrating, in an adaptive filter, parameters relating to lateral dynamics and 
10 parameters relating to longitudinal dynamics of said vehicle; 

estimating the values of said integrated parameters dependent on a lateral 
movement indication signal and on a road friction indication signal derived from a 
wheel angular velocity signal, a motor torque signal and a gear indicator signal; 
calculating an improved tire pressure indication signal for an identified tire, and 
15 possibly also a confidence indication, dependent on said integrated parameter values. 

7. The method as recited in any of the preceding claims, wherein said lateral movement 
indication signal alternatively or additionally is generated by a sensor indicating the 
lateral acceleration of said vehicle, said sensor preferably being an accelerometer. 

20 

8. The method as recited in any of the preceding claims, further being adapted for tire 
pressure estimation in a motorcycle having two wheels by comprising in said adaptive 
filter front and rear wheel angiilar velocity parameters. 

25 9. The method as recited in any of the preceding claims, further being adapted for tue 

pressxire estimation in a motorcycle having two wheels by comprising in said adaptive 
filter a roll angle parameter, the value of which is represented by a roll angle sensor 
signal. 

30 1 0. The method as recited in claim 1, wherein said adaptive filter is based on a 

predetermined model of vibrations in a tire, and said first tire pressme indication value 
is a resonance frequency. 

1 1 . The method as recited in claim 1 , wherein the adaptive filter of said predetermined 
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model is a recursive least squares algorithm. 

12. The method as recited in claim 1, wherein the adaptive filter of S2iid predetermmed 
model is a Kalman filter representing a state space model algorithm. 

5 

13. The method as recited in claim 1, further comprising the step of pre-processing a 
vehicle status signal in the shape of a wheel angular velocity signal derived from a 
toothed wheel sensor with respect to correction of toothed-wheel errors. 

10 14. The method as recited in claim 13, further comprising the step of interpolating non- 
equidistant data in said corrected wheel angular velocity signal. 

15. The method as recited in any of the preceding claims, fiu1:her comprising the step of 
reducing the sampling frequency of a sampled sensor signal by re-sampling and band 

15 pass filtering the input sensor signals into a first and a second frequency band, and by 
modulating said first and second frequency bands with modulating signals dependent 
on ttie lower frequency of each frequency band. 

16. The method as recited in claim 1, fiirther comprising the step of: 

20 integrating, in an adaptive filter, parameters relating to different groups of vehicle 

status parameters; 

estimating the values of said integrated parameters dependent on a sensor signals 
representing entities related to said groups of vehicle status parameters; 

calculating a tire pressure indication signal dependent on said integrated parameter 
25 values, and possibly also cedculating secondary indication signals derived from the 
values of said integrated parameters. 

17. The method as recited in claim 1, further comprising the step of: 

calculating a first estimate of a tire air pressure value dependent on said first tire 
30 pressure indication value and a predetermined mapping function or mapping table. 

18. A computer program product for estimating the pneumatic pressure of a tire on a 
veliicle, comprising program code adapted to direct a data processing system to 
perform the steps and functions of the preceding claims. 
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19. An apparatus for estimating the pneumatic pressure of a tire on a vehicle, comprising 
functional components for perforaiing the steps and functions of any of the preceding 
claims. 
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